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We present a method for achieving complete population transfer in a two-state quantum system 
via adiabatic time evolution in which, contrary to conventional rapid adiabatic passage produced 
by chirped pulses, there occurs no crossing of diabatic energy curves: there is no sign change of 
the detuning. Instead, we use structured pulses, in which, in addition to satisfying conditions for 
adiabatic evolution, there occurs a sign change of the Rabi frequency when the detuning is zero. We 
present simulations that offer simple geometrical interpretation of the two-dimensional motion of 
the Bloch vector for this system, illustrating how both complete population inversion and complete 
population return occur for different choices of structured pulses. 

PACS numbers: 32.80.Xx, 33.80.Be, 32.80.Qk, 03.65.Ge 



I. INTRODUCTION 



(RWA) [J, (2 reads 



The technique of rapid adiabatic passage (RAP) [l|, 
& 1 S S provides a well-studied robust and efficient 
method for producing complete population transfer be- 
tween two bound states of a quantum system. Tradi- 
tionally the population change has been accomplished by 
sweeping the carrier frequency of a laser pulse through 
resonance with a two-state system (parameterized by 
detuning), while simultaneously pulsing the strength of 
the interaction (parameterized as time-varying Rabi fre- 
quency) . The sweep of detuning corresponds to a cross- 
ing of diabatic energy curves, as noted below. When the 
interaction changes sufficiently slowly, so that the time 
evolution is adiabatic, there will occur a complete popula- 
tion transfer from the ground state (of energy E±) to the 
excited state (energy E2). Alternatively, one can make 
the detuning cross resonance by sweeping the transition 
frequency, e.g. by time-dependent electric or magnetic 
fields, as exemplified by the recently introduced tech- 
nique of Stark-chirped rapid adiabatic passage (SCRAP) 
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In this paper we exploit a novel way of achieving com- 
plete population transfer via adiabatic evolution in two- 
state systems, using pulses that produce no crossing of 
diabatic energies. Such population transfer is possible 
when diabatic energy curves touch, but do not cross, and 
simultaneously there occurs a sign change of the interac- 
tion (the Rabi frequency). 



II. BACKGROUND 

Coherent excitation of a two-state quantum system 
is described by the time-dependent Schrodinger equa- 
tion (TDSE), which in the rotating- wave approximation 



m-c(t) = H(«)c(t), 



(i) 



where C(t) 



[Ci(t),C2(t)] T is a two-dimensional 



column- vector whose elements are the probability ampli- 
tudes Ci(t) and C2{t) associated with the two quantum 
states ipi and ip 2 , and H (t) is the matrix representing the 
RWA Hamiltonian 0,0], 



-A(t) 0(t) 

n(t) A{t) 



(2) 



The detuning A(i) = Lu (t) — wl(£) expresses the fre- 
quency offset of the laser carrier frequency wl(£) from 
the Bohr transition frequency wo(i), each of which may 
vary with time. The Rabi frequency f2(t) quantifies the 
field-induced coupling between the two states. For the 
laser-atom excitation considered here this is the product 
of the atomic dipole transition moment d\2 in the field 
polarization direction, and the electric-field envelope £(t) 
of the laser field, fl(t) = —di^Eit) jh. We shall take this 
to be a real-valued function of time. 

In writing the RWA Hamiltonian we have chosen a 
symmetrized form, which is obtained by a suitable phase 
transformation of the probability amplitudes. The di- 
agonal elements of H(t) are diabatic energies. Curves 
of diabatic energies cross when A(t) changes sign. The 
(time-dependent) eigenvalues of H(t) arc adiabatic ener- 
gies, He±(t). By definition, £+(£) is the larger of the two 
eigenvalues. 

As has long been known 0, Eq ([lj) for the two 
complex-valued probability amplitudes can be recast as 
three coupled equations for real- valued variables, that 
serve as the three components of a unit vector B(i) = 
[u(t),v(t),w(t)] T in a three-dimensional abstract space. 
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The components of this Block vector are 

u(t) =2Re p 12 (t), 
v(t) = 2Impi 2 (f), 

W(t) = P22{t) - pil{t), 



(3a) 
(3b) 
(3c) 



with p mn (t) = C* n (t)C n (t) (n, m = 1,2) being the matrix 
elements of the density operator. The resulting RWA 
optical Bloch equation reads [l|, 0, 0] 



— B(t) = R(t)B(t), 



with 



R(t) = 





AW 




-A(i) 




o 

-f2(t) 




(4) 



(5) 



The skew-symmetric nature of the matrix R(t) allows us 
to write Eq (j4|) as a torque equation for the Bloch vector 

0,9, 



at 



B(t) = Q(t) x B(t), 



(6) 



where the driving torque (angular velocity) vector Q(i) 
reads @,H 

Q(i) = [r!(i),0,A(t)] T . (7) 
The Bloch vector rotates at the instantaneous rate 



n(t) = v/f>(t) 2 + AW 2 . 



(8) 



Note that the direction of rotation will be reversed if the 
signs of f2(t) and A(t) are both reversed. 

One of the three eigenvectors of R(t), which corre- 
sponds to a null eigenvalue, is a sum of only two compo- 
nents, the coherence u(t) and the inversion w(t), 



a (t) = u(t) sintf(i) +w(t) costf(i). 



(9) 



It has no component of v(t); in analogy with stimulated 
Raman adiabatic passage (STIRAP) [1, 0, [HI, we shall 
refer to it as a dark vector. Here the time-dependent 
mixing angle i)(t) is defined by 



tantf(t) = n{t)/A(t). 



(10) 



In contrast to the dark state of STIRAP, the dark vec- 
tor ao(t) can radiate via spontaneous emission from the 
excited state ip2 ■ It is important to note that the compo- 
nents of the dark vector are the normalized components 
of the torque vector ([7]) . 

Below we describe several situations in which the time 
evolution is adiabatic, so that the Bloch vector adiabat- 
ically follows the torque vector as, in accord with the 
design of the functions A(t) and 0(t), it moves in a 
two-dimensional plane of the abstract vector space. The 
pulses can be designed to produce either complete popu- 
lation inversion (CPI), in which all population resides in 
state 2 at the completion of the pulses, or complete pop- 
ulation return (CPR), in which the initially populated 
state is completely rcpopulated after the pulses. 
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FIG. 1: CPI (left frames) and CPR (right frames) in the 
presence of level crossing (detuning crossing resonance) . Left 
frames: The detuning and the Rabi frequency are given by 
Eqs CED, with A = -40 /T and fi = 30/T. Right frames: 
the detuning and the Rabi frequency are given by Eqs (|12|l . 
with r = 3T, A = -25/T and Q = -80/T. Top frames: 
the Rabi frequency and the detuning; second frames: eigenen- 
ergies; third frames: mixing angle $(t); bottom frames: com- 
ponents of the Bloch vector B(t). 



III. ADIABATIC FOLLOWING WITH LEVEL 
CROSSING 

Figures [1] shows examples of adiabatic CPI or CPR in 
the presence of a level crossing (the detuning crosses the 
resonance) . The frames in the left column of Fig. [T] il- 
lustrate CPI for the conventional rapid adiabatic passage 
(RAP) by a frequency chirped laser pulse [H, S, S, 0, Si • 
Top frames show the pulsed Rabi frequency and detuning 
as a function of time, for which we have chosen 



A(i)^e-W L5T > ! 

n(t) = n e-^ 2 , 



(11a) 
(lib) 



The second frames show the eigenencrgies, the third 
frames show the mixing angle $(t), and the bottom 
frames show the components of the Bloch vector. 

The frames in the right column of Fig. Q] illustrate the 
situation when, in addition to the level crossing of the 
detuning, the Rabi frequency also changes sign, 



A(t) = 



Apt 2 
T 2 



sech 



(t/T) 2 



t ~ T 

T 



— sech 



t + T 

T 



1,2a) 
(12b) 



The temporal pulse area - the time integral of the Rabi 
frequency - is here zero (a "zero-area pulse"). For this 



3 




-10 -5 5 
Scaled Time t/T 



10 -5 5 10 
Scaled Time t/T 



FIG. 2: CPI (left frames) and CPR (right frames) in the 
absence of level crossing (detuning only touching resonance). 
Left frames: the detuning and the Rabi frequency are given 
by Eqs with r = 3T, A = 4/T and fl = -40/T. 

Right frames: the detuning and the Rabi frequency are given 
by Eqs JUJ, with r = 3T, A = 4/T and ft = 40/T. Top 
frames: the Rabi frequency and the detuning; second frames: 
eigenenergies; third frames: mixing angle $(t); bottom frames: 
components of the Bloch vector B(t). 



IV. ADIABATIC FOLLOWING WITHOUT 
LEVEL CROSSING 

Figure [5] shows a pair of cases, which do not appear 
to have been reported previously. In the left frames of 
Fig. [2] we illustrate CPI in the absence of a level crossing, 
when the detuning does not cross but just touches the 
resonance. The detuning and the Rabi frequency read 



A(i) = 



sech 



fi(t) = -£-e 



Apt 2 

J^2 



-(t/2.5T) 2 



T 



sech 



t - T 
T 



1,3a) 
(13b) 



During the excitation the mixing angle i}(t) rotates from 
i9(— oo) = to #(oo) = 7T, which leads to CPI in the adia- 
batic limit because the system remains in the dark vector 
Furthermore for the given example of antisymmet- 
ric pulse shape, the time-integrated Rabi frequency (the 
temporal pulse area) is zero, but the zero pulse area is 
not a necessary constraint. This is a new approach to 
CPI, without a level crossing and with a possible zero 
pulse area. 

The right frames of Fig. [5] illustrate the case when the 
Rabi frequency and the detuning are both even functions 
of time, 



fi(0 = tt e- (t/2 - 5T): 



t + T 

T 



sech 



t - T 
T 



1,4a) 
(14b) 



choice of antisymmetric Rabi frequency and detuning, 
exact CPR takes place. Indeed, this property, which 
stems from the overall antisymmetry of the Hamiltonian, 
is valid not only in the adiabatic limit but in the general 
case, and for arbitrarily many states [l2T |. 

For the left frames, during the excitation the mixing 
angle i9(t) rotates from i9(— oo) = to $(0) = tt/2 and 
then to -d(oo) = ir (due to the change of the sign of the 
detuning), which leads to CPI because the system re- 
mains in the dark vector ©, which follows adiabatically 
the driving torque. For the right frames, during the ex- 
citation the mixing angle $(t) rotates from •&(— oo) = 
to i?(0) = 7r/2 and then back to i?(oo) = (due to the 
change of sign of the detuning and the Rabi frequency at 
the crossings), which leads to CPR because, again, the 
system remains in the dark vector (O, but now it has 
different asymptotics. 

The choice of peak amplitudes for the pulses in these 
examples is such that a measurable value of Bloch vector 
component v(t) occurs during intermediate stages of the 
process. As the motion becomes more adiabatic, and 
the adiabatic following more complete, this component 
becomes smaller. 



The mixing angle #(t) rotates from oo) = to $(0) = 
7r/2 and then to i?(oo) = 0, which leads to CPI in the 
adiabatic limit because the system remains in the dark 
vector ©. 



V. DISCUSSION 
A. General features 

In summary, the analysis above brings us to the follow- 
ing conclusions. CPI takes place in the adiabatic limit 
when (i) the coupling is an even function of time and 
the detuning is odd function of time (f2(t) = fi(— t), 
A(t) = —A(—t), left frames of Fig. [J), or (ii) the coupling 
is odd function of time and the detuning is even function 
of time (n(i) = -fi(-i), A(t) = A(-i)), left frames of 
Fig. [2 When (iii) both functions are even function of 
time or (iv) both are odd function of time, then CPR is 
observed instead of CPI, right frames of Figure [T] and [2] 
These findings are summarized in Table [T] 

We emphasize here that the assumption of definite 
symmetry of the Rabi frequency f2(t) and the detuning 
A(t), which was made in figures [T] and [21 and summa- 
rized in Table [H is only a matter of convenience. In the 
adiabatic following approximation, which was assumed 
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Figure 


Rabi frequency 


detuning 


adiabatic 




n(-t) 


A(-i) 


limit 


U (left) 




-A(i) 


CPI 


[H (right) 


-n(t) 


-A(t) 


CPR 


M (left) 


-n(t) 


AW 


CPI 


[3 (right) 




A(f) 


CPR 



TABLE I: Adiabatic limit for fields of particular symmetry: 
CPI: complete population inversion; CPR: complete popula- 
tion return. 




throughout, it is only the asymptotic behavior of the mix- 
ing angle $(t), which controls the composition of the dark 
vector that matters. 

Because the adiabatic passage is robust, the procedure 
of CPI with level touching (left frames of Figure [2]) is 
also robust: it depends only weakly on the overlap of the 
coupling fi(f) and the detuning A(t) and their peak val- 
ues. Moreover, the exact antisymmetry is not required, 
and high efficiencies can also be achieved with asymmet- 
ric time dependences. The only restriction of the pattern 
that we have to follow in order to achieve CPI is: the 
detuning starts before the coupling, becomes zero before 
the coupling changes sign, then rises again with the same 
sign and ends after the coupling ends. 

The torque-equation interpretation of the motion of 
the Bloch vector provides a useful means of understand- 
ing the behaviors discussed above. Figure[3]provides such 
a picture for the four frames of figures Q] and [2] and Ta- 
bled] In the adiabatic limit, which we assume, the Bloch 
vector evolves in the plane (it, w), which is the v — slice 
of the Bloch sphere. The geometrical representation of 
the adiabatic evolution of the Bloch vector is easily de- 
termined because in the adiabatic limit the state vector 
coincides with the dark vector ([9]) and is always parallel 
or anti-parallel to the driving torque. 



B. Experimental feasibility 

The concept of a zero-area laser pulse, when first en- 
countered, may seem somewhat confusing. We have dis- 
cussed such pulses in some detail elsewhere [IH . Here we 
describe briefly a few of the available techniques. 



1. Nonoscillating magnetic fields. 

Historically, the concept arose first in nuclear magnetic 
resonance (NMR), wherein typically spin-^ particles are 
controlled by slowly varying magnetic fields fl4j . A zero 
area field can be produced in an obvious manner, by ap- 
propriate tailoring of an applied voltage that can reverse 
the direction of the field and create the desired asymmet- 
ric pulse. In NMR the detuning A(t) is replaced by a lon- 
gitudinal (along the quantization axis) magnetic pulse, 
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FIG. 3: Motion of the Bloch vector and the angular veloc- 
ity vector in the (u, w) plane. Solid line (blue online) is the 
angular velocity (torque) vector, normalized to peak value. 
Dashed line (red online) is the Bloch vector. The initial posi- 
tions of the Bloch vector and the torque are marked by hollow 
circles. The final positions are marked by filled circles. The 
cases (a)-(d) correspond to the cases of figures [T] and [2] 

and the Rabi frequency il(t) by a transverse magnetic 
pulse. 

2. Self-induced transparency. 

In the optical spectrum zero-area pulses arc formed 
when pulses are reshaped by self-induced transparency 
(SIT) while propagating through an absorbing medium. 
Initial pulse areas in the range (0, 2ir) will tend toward 
zero-area pulses - an expression of the McCall-Hahn area 
theorem |l5|. 

3. Beam splitting & recombination. 

A zero-area pulse can also be produced by passing a 
laser pulse through a beam splitter and then recombining 
the two split pulses, after delaying one of them by a time 
delay r and shifting its phase by it. The recombincd 
pulse, with Rabi frequency 

n(t) = n [f(t-r)-f(t)}, (15) 

where r is the delay, is a zero-area pulse. Such an im- 
plementation requires an interferometric stability within 
a wavelength A. 
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4- Ultrashort shaped pulses. 

The most useful techniques, for ultrashort (femtosec- 
ond and picosecond) laser pulses, manipulate the Fourier 
transform of the pulse, 

~ i r°° 

ft(A) = — dt e iAt fl(t), (16) 

using pulse shapers [lj|. The pulse area A is the value 
of this transform on resonance (A = 0), 

/oo 
dtn(t) = 2irtt{0). (17) 
-oo 

Thus a zero-area pulse can be produced by negating the 
resonant frequency component 0(0) in its spatially dis- 
persed Fourier spectrum. 

The desired time variation of the detuning can be pro- 
duced, using a pulse shaper, simultaneously with the ma- 
nipulation of the Rabi frequency. Alternatively, the pulse 
of Fig. [^c) can be produced by shaping the Rabi fre- 
quency and the detuning separately. For example, the 
Rabi frequency can still be shaped by a pulse shaper, 
whereas the detuning can be shaped by manipulating the 
Bohr transition frequency (rather than the laser carrier 
frequency through the spectral phase) by using a pair 
of suitably delayed far-off-resonant Stark-shifting laser 
pulses, as in SCRAP [||. 

C. Comparison with earlier work 

CPI with a zero-area pulse has been reported in previ- 
ous publications [17| . However, the physical mechanisms 
in these publications is different from the one reported 
here because the CPI in the present paper relies on adia- 
batic time evolution at all times. By contrast, in the pre- 
vious proposals [l7| the sudden change of the phase of the 



field was essential, which led to sudden evolution, in the 
form of a delta-function shaped nonadiabatic coupling in 
the adiabatic basis. This required the use of pulses with 
large Rabi frequencies, much larger than here. 

It is important to note that the zero-area condition 
is not required in any of these papers because the CPI 
mechanisms do not require it. The completely antisym- 
metric shape of the Rabi frequency, f2(— t) = — ft(t), and 
the ensuing zero pulse area, is the most natural realiza- 
tion of a function that changes sign at time t = 0, but it 
is by no means necessary for the physical mechanisms. 



VI. CONCLUSIONS 

We have shown that complete population transfer, by 
means of adiabatic passage, is possible without any cross- 
ing of diabatic energy curves. It is necessary only that 
the curves touch, and that there be a sign change in the 
Rabi frequency at that moment. Such a scenario could be 
used to replace the conventional two-state rapid adiabatic 
passage via curve crossing for various implementations. 
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